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ABSTRACT 


Fallout  mass-size  distributions  presently  used  at  USNRDL  are  com¬ 
pared  with  new  distributions  suggested  by  recent  investigations.  Avail¬ 
able  data  is  unable  to  define  the  distribution  parameters  well  enough 
to  distinguish  between  lognormal  and  power-lav  distribution  models. , 


SUMMARY 


Problem 


To  determine  the  differences  between  newly-proposed  fallout  mass- 
size  distributions  and  distributions  now  in  use. 

Findings 

Differences  are  trivial  in  comparison  with  the  effects  of  current 
uncertainties  in  the  distribution  parameters.  Planned  sensitivity 
analyses  need  not  explicitly  treat  the  newly -proposed  power-law  distri 
buttons  except  as  a  matter  of  convenience. 
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INTRODUCTION 


The  purpose  of  this  report  is  to  compere  formulas  for  the  size 
distribution  of  particle  mass  in  the  fallout  from  land-surface  bursts. 

The  distribution  equations  to  be  compared  are  basically  the  newly- 
proposed  power-law  distribution  and  the  lognormal  distribution,  although 
modifications  of  each  are  involved.  Cumulative  forms  of  the  distribu¬ 
tions  will  be  emphasized.  Interest  in  such  a  comparison  arises  from  the 
need  to  assess  the  impact  of  new  information  and  suggestions  (mainly  due 
to  Russell)1  on  present  prediction  techniques.  The  mass-size  distribu¬ 
tion  is  of  basic  importance  for  predicting  fractionation  effects  because, 
being  equivalent  to  the  volume -size  distribution,  it  is  primary  input 
data. 


At  this  Laboratory  both  total  particle  mass  and  active  particle 
mass  are  considered  to  have  the  same  lognormal  distribution,  even  though 
there  may  be  some  fifteen  times  as  much  inactive  material  as  active 
material. 2  The  power-law  distribution,  proposed  by  Russell,1  refers 
only  to  total  mass.  The  land-surfaces  involved  are  both  silicate 
(Nevada  Test  Site -NTS)  and  coral  (Eniwetok  Proving  Grounds -EPG) . 

The  results  of  the  mathematical  comparison  offered  here  will  throw 
light  upon  the  need  of  a  sensitivity  analysis  for  more  detailed  compari¬ 
son. 
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BACKGROUND 


For  purposes  of  fallout  prediction,  the  information  required  for 
input  can  be  conveniently  discussed  in  terms  of  cumulative  logarithmic- 
probability  graphs  because  any  single -valued  distribution  curve  on  such 
a  graph  is  automatically  normalized  to  100  Let  us  say  that  the  pre¬ 
dictor  is  interested  in  the  distribution  of  Mo99  in  a  local  fallout 
field.  Even  though  his  computer  can  reasonably  handle  hundreds  of  par¬ 
ticle-size  classes,  the  data  on  vhich  his  input  is  based  need  not  be 
so  detailed.  Presumably  he  knows  enough  about  the  device  to  state  the 
total  quantity  of  Mo99  present.  His  next  most  urgent  need  is  to  know 
vhat  fraction  of  this  came  down  locally,  or  at  early  times,  or  in  large 
particles.  This  corresponds  to  placing  a  point  on  the  graph  in  the 
region  of  25  to  50 -n  diameter  particles.  Such  a  point  might  conceivably 
come  from  several  sources,  e.g., 

1.  Prediction  techniques  in  vogue, 

2.  Partition  inferred  from  the  radiochemical  analysis  of  cloud 

and  ground  samples,  gq 

3*  Integration  of  contours  of  Moy^  surface  density  on  the  ground. 

Assuming  a  number  is  available  to  place  this  point,  he  still  needs 
199  more  numbers.  These  can  be  obtained  from 

1.  Prediction  techniques  in  vogue, 

2.  Determinining  Mo99  distribution  with  particle  size  in  a  repre¬ 
sentative  sample  of  the  local  field, 

3*  Integration  of  contours  of  particle  surface  density  as  a  func¬ 
tion  of  size,  plus  a  knowledge  of  the  relative  Mo99  content  of  different 
particle -size  classes, 

k.  Determining  Mo99  distribution  with  particle  size  in  a  cloud 
sample  which  contains  representative  proportions  both  of  the  particles 
chosen  above  and  of  larger  particles. 

The  more  observational  data  that  can  be  incorporated,  the  less  reliance 
need  be  placed  on  choice  one  and  moreover,  the  better  the  quality  of 
choice  one  that  can  be  offered.  If  as  few  as  five  well -distributed 
points  were  available  from  observations,  the  predictor  would  be  in  good 
shape. 
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These  same  considerations  apply  to  the  distribution  of  other  par¬ 
ticle  properties:  individual  radionuclides,  mass,  or  fraction  of  the 
unpartitioned,  unfractionated  normalization  factor.  Unfortunately,  it 
is  invariably  necessary  to  rely  on  choice  one  to  some  extent.  It  is 
therefore  incumbent  upon  the  predictor  to  keep  current  the  comparison 
of  his  techniques  with  published  observations. 


THE  LOGNORMAL  DISTRIBUTION 


To  illustrate  the  notation  we  may  introduce  the  lognormal  distribu 
tion  by  saying  that,  if  the  probability  is  (a,b)  that  a  randomly 
chosen  particle  in  a  distribution  has  a  mass  x:  a  <  x  <  b,  then  the 
equations  ^ 

PM  <a’tJ  ■  f  %  M  dx  I 

a 

PM  {x>x+dxi  =  %  M  dx  J 

define  the  probability  density  function  p„(x),  and  for  the  lognormal 
distribution. 


exp  | 

The  parameters  u,  and  c  are  best  determined  from  a  cumulative  log-prob¬ 
ability  plot  of  mhe  mass  distribution.  If  a  straight  line  gives  a 
reasonable  fit  to  the  cumulative  weight,  then 

Hh  =  ln 

(where  x^  is  the  particle  diameter  at  50  $) 


0  =  In  \+w  ln  \ 

V 
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where  x^+  and  x^-  are  the  particle  diameters  one  sigma  unit  larger  and 
smaller  than  34  (i.e.,  at  84.13  $  and  15*87  %  respectively).  The  mass 
distribution  currently  in  use  at  NRDL  uses  the  lognormal  distribution 
with  parameters  xjj  =  100  4  and  a  «  1.68?. 

The  mas s-  distribution  curve  favored  at  NRDL  is  based  upon  (1)  the 
assumption  of  constant  volume -specific  activity  and  (2)  activity-size 
data  from  75  to  3300  4. 2  An  argument  can  be  made  for  assuming  constant 
surface-specific  activity  instead  of  constant  volume- specific  activity. 
The  effect  of  this  will  be  discussed  below. 


THE  POWER  DISTRIBUTION 


In  an  unclassified  section  of  a  classified  report,*  Russell  des¬ 
cribes  his  analysis  of  Johny  Boy  cloud  samples,  which  analyses  led  him 
to  the  conclusion  that  the  mass  of  the  debris  was  so  distributed  that, 
down  to  about  90  4,  equal  size  increments  contained  equal  masses: 

PM  (a,b)  «(b-a) 


or 


=  constant 


On  the  assumption  of  constant  particle  density  and  the  usual  spherical- 
particle  approximation,  this  leads  to  a  frequency  distribution  function 

PN  {x,x-hix}  ==  Pjj  (x)dx  =  lyc^dx 

where  kjj  is  a  normalization  factor. 

Further  investigation,  especially  by  Nathans,*  indicated  variations  from 
the  power  of  three,  so  for  generality  ve  write 

PN  {x,x-KJx)  =  kjjx”%x 

PM  {x,x+dx}  =  ^(x)  dx  =  k^x3"^  dx 

which  bears  little  resemblance  to  Equation  1.  The  value  of  q  appears 
to  lie  between  3  and  4. 


*M.  W.  Nathans,  Tracerlab,  Inc.,  private  communication. 
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Completion  of  the  distribution  equation  requires  that  some  upper 
limit  Xjjjqx  be  set  to  the  particle  size.  Hius  the  form  of  the  equation 
becomes  ,  , 

b  ,  „  /jw  ,  n 

PM  ia'bi  “  J  9  dx/ J  x3  q  dx 

a  /  o 

Various  considerations  led  Russell  to  choose  a  value  of  1000  p  for  xMX, 
although  References  2  and  3  indicate  a  higher  value  would  be  more 
appropriate . 

Figure  1  compares  the  power  distribution  with  the  lognormal  distri¬ 
bution  for  several  cases.  The  power  distribution  is  plotted  for  typical 
values  of  q  =  3*0  and  3»5>  both  with  a  cutoff  of  x^y  =  1000  p.  A 
curve  for  q  =  3-0  and  =  2000  p  shows  the  effect  of  a  reasonable 

change  in  cutoff  diameter. 

For  the  case  of  q  =  3*0  and  xmHY  =  1000  p  the  mass  distribution  is 
extremely  simple.  The  fraction  of  mass  in  particles  of  diameter  x  or 
less  is  then 

p  /o.x}  « 

M  *  1000 


For  the  case  of  q  =  3*5  and  x  =  1000  p  it  is  not  much  more  com- 
,  .  .  .  max 

plicated : 

x  / 1000 


PM  {0,x>  =  |  x"l/2  dx  / f  x"l/2 


1000 


For  the  lognormal  distributions  the  equations  are  derived  according 
to  the  method  of  Reference  4 : 


A 


■n?.)23-°P.3  d  In  x 

•  6o2  j 


for  constant  volume -specific  activity  and 

x  / 


In  x 


for  constant  surface-specific  activity. 
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PARTICLE  DIAMETER,  x(M) 


PERCENT  OF  MASS  IN  PARTICLES  SMALLER  THAN  SIZE  INDICATED 


Fig.  1  Comparison  of  Extended  Distributions 
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THE  TRUNCATED  DISTRIBUTIONS 


The  Truncated  Power  Law 


It  was  immediately  evident  to  Russell  that  the  power -law  distribu¬ 
tion  was  physically  unrealistic  as  shown.  Thus  as  x  -*  0,  the  distribu¬ 
tion  predicts  an  infinite  number  of  particles.  Therefore  Russell  chose 
a  value  of  1  p  for  the  lower  cutoff  diameter  x^n,  giving  the  final 
form  of  the  distribution  as  b  A"! 

f 

}1000  dx 

1 

The  small  size  cutoff  has  since  been  well  documented  by  Nathans’ 

vork . 

The  Truncated  Lornormal  Distribution 


If  pressed  to  the  limit,  the  lognormal  distribution  will  predict 
some  finite  (if  fractional)  number  of  particles  of  any  size,  no  matter 
how  large  or  small.  A  means  of  truncating  the  lognormal  distribution 
at  the  small  particle  sizes  x  .  was  employed  in  Reference  4.  This 
involved  simply  making  the  substitution 


x-x 


'min 


in  the  distribution  function.  If  a  cutoff  at  larger  sizes  is  desirable, 
one  needs  another  means  of  making  the  logarithmic  term  infinite.  This 
can  be  easily  supplied  by  first  writing 


In  x  -  -  In  - 
x 


and  making  the  substitution 


max 


which  is  equivalent  to  making  the  substitution 


max 


I 


I 


An  obvious  way  to  incorporate  both  cutoffs  simultaneously  is  to  make 
the  single  substitution 


X  -  X 


min 


1  -  x/j 


max 


Appendix  A  describes  how  truncated  graphs  can  be  constructed  from  the 
parameters  and  o  with  a  minimum  of  effort. 

Comparison 


Figure  2  compares  the  lognormal  and  power -lav  distributions  in 
truncated  form.  For  comparison,  all  distributions  are  truncated  at  1  p 
and  1000  p.  With  the  truncated  power-law  distribution  it  is  possible 
to  show  the  graph  for  q  =  4.  Therefore,  Fig.  2  shows  the  full  extent 
of  uncertainty  due  to  the  uncertainty  in  q.  FOr  the  case  of  q  =  4: 


d  In  x 


It  is  noteworthy  that  this  q  =  4  distribution  is  also  the  distribution 
of  surface  that  would  correspond  to  the  q  =  3  mass  or  volume  distribution. 


The  truncated  lognormal  distributions  are  shown  for  the  assumptions 
of  both  constant  volume-specific  activity  and  constant  surface-specific 
activity.  The  former  lies  well  within  the  range  of  power-law  distribu¬ 
tions  . 


Variation  in  x 

max 


As  indicated  above,  the  value  of  Xjnax  fs  not  well  known.  The  above 
equations  are  easily  adapted  to  other  values.  Choosing  500  and  2000  for 
illustrations  only  gives  the  values  shown  in  Table  1.  The  sensitivity 
to  Xmfly  is  seen  to  be  strongly  dependent  on  the  value  of  q.  Thus,  the 
effect  on  %  {0,x}  ranges  from  a  factor  of  1.2  for  q  =  4  up  to  a  factor 
of  four  for  q  =  3* 


8 


NMOL  71*-** 


PERCENT  OF  MASS  IN  PARTICLES  SMALLER  THAN  SIZE  INDICATED 


Fig.  2  Comparison  of  Truncated  Distributions;  Effect  of  Uncertainty  in 
Exponent 
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Effect  of  Varying  xmax  on  Power -Law  Equations  with  Various  Exponents 


Cutoff  Diameter 

PM 

<0,x> 

q  -  3-0 

q  = 

q  «  4.0 

500  u 

x-l 

V™. 

1  , 

E99 

•■iir 

2.7  logio 

1000  u 

X'l 

1  , 

3o75~ 

3  logiox 

2000  \i 

X-l 

V^i 

1  1 

1999 

53T T 

3.3  logio  x 

SIGNIFICANCE 


As  Fig.  1  illustrates,  for  the  particle  size  range  of  10-800  \i,  the 
extended  power-law  distributions  for  exponents  of  3*0  and  3*5  lie  between 
the  lognormal  curves  based  upon  constant  surface-specific  activity  and 
constant  volume -specific  activity.  Below  10  p,  the  divergence  of  the 
power-law  distribution  from  the  area  bounded  by  the  lognormal  lines  is 
at  most  about  k  $  of  the  total  mass.  Above  800  \x,  this  divergence 
varies  with  large  particle  cutoff,  and  for  a  cutoff  of  2000  ji  the  diver¬ 
gence  is  less  than  5  #. 

Figure  2  shows  that,  for  truncated  distributions,  the  lognormal 
curve  based  on  the  assumption  of  constant  volume -specific  activity  lies 
well  within  the  range  of  uncertainty  in  power-law  exponents  and  does  not 
differ  greatly  fran  the  curve  based  on  an  exponent  of  3»5-  Although  the 
lognormal  curve  for  the  assumption  of  constant  surface -specific  activity 
lies  outside  the  range  of  uncertainty  in  power-law  exponents,  it  does 
not  differ  greatly  from  the  curve  based  on  an  exponent  of  3‘0- 

Figure  3  is  a  striking  illustration  of  the  similarity  that  can  be 
achieved  between  the  truncated  power-law  distribution  and  the  truncated 
lognormal  distribution.  This  figure  compares  the  q  =  3.5  power-law 
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PARTiCLE  O  AMETER,  x(^t) 


) 


NROL  716  66 


PERCENT  OF  MASS  IN  PARTICLES  SMALLER  THAN  SIZE  INDICATED 


Pig.  3  Achievable  Similarity  Between  Power-Law  and  NKDL  Lognormal 
Distributions.  Lognormal  parameters  are  =  500  \x,  x^.  =  30  p, 

xmin  '  1  ^  W  1000  •*' 
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distribution  with  the  ^  =  500  n,  x^j-  =  jO  p  lognormal  distribution. 

The  agreement  is  such  that  if  one  curve  is  correct,  the  other  will 
never  be  proved  wrong.  Presumably,  equal  similarity  could  be  achieved 
with  curves  for  other  values  for  q.  Presumably  also,  the  agreement  in 
Pig.  3  could  be  made  still  closer  by  either  (a)  a  different  choice  of 
a  and  |^|  or  (b)  an  equally  reasonable  method  of  cutting  off  the  log¬ 
normal  distribution  at  1000  p,  or  (c)  cutting  off  the  lognormal  distri¬ 
bution  at  some  equally  defensible  value  like  950  p. 

Thus  it  appears  that  th_*  differences  between  the  two  approaches 
are  trivial.  The  lognormal  distribution  has  the  esthetic  advantage  of 
an  observationally  confirmed  theoretical  basis  in  the  case  of  airburst 
debris.  If  truncation  is  required,  the  power-law  distribution  has  the 
practical  advantage  of  simplifying  further  calculations  of  particle 
surface  distribution.  Planned  sensitivity  analyses  need  not  incorpor¬ 
ate  power-law  distributions  explicitly,  but  may  include  them  as  a  matter 
of  convenience. 
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APHND3X  A 


GRAPHICAL  CONSTRUCTION  OP  TRUNCATED  LOGNORMAL  DISTRIBUTIONS 


Consider  first  the  region  of  sizes  less  then  x  and  assume  that  the 
effect  of  Xj,^  is  negligible  in  this  region.  The  truncated  curve  can 
be  calculated  from  the  extended  curve  by  simply  replotting  percentages. 
Thus,  for  a  cutoff  of  1  p,  the  value  of  the  extended  curve  for  1  p  is 
replotted  at  2  p,  the  value  for  the  extended  curve  at  2  p  is  replotted 
at  3  p,  etc. 

In  the  size  region  above  x,  assuming  the  effect  of  x_£_  can  be  neg¬ 
lected,  the  procedure  is  similar.  Thus,  for  a  value  of  1000  p  for  x^^., 
if  one  wishes  to  plot  a  value  for  900  p,  one  simply  plots  the  value 
corresponding  to  900/(1  -  900/1000)  or  9000  p  on  the  extended  curve. 
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